General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



FACILITY FORM 602 


ACOUSTIC SCATTERING BY A POROUS ELLIPTIC CYLINDER 


WITH NONLINEAR RESISTANCE 


By 


William Edward Zor\jmski 


Thesis submitted to the Graduate Faculty of the 
Virginia Polytechnic Institute 

in candidacy for the degree of 

DOCTOR OF PHILOSOPHY 


in 




Engineering Mechanics 

8 ^ 


" 3 ^ 


(CODE) 


(CATEGORY) 




March 1970 



ACOUSTIC SCATTERING BY A POROUS ELLIPTIC CYLINDER 


WITH NONLINEAR RESISTANCE 

by 

William Edward Zorxir"ki 

•Hiesis submitted to the Graduate Faculty of the 
Virginia Polytechnic Institute 
in partial fulfillment for the degree of 

DOCTOR OF PHTLOSOPHY 

in 

Engineering Mechanics 


APPROVED: 



Dr. ft. 'rtiomas Davis 











Dr. Dean T. Mook 


Dr. George A. Gray 




iy 


March 1970 


Blacksburg, Virginia 


ACOUSTIC SCATTERING BY A FOROUS ELLIPTIC C^fLINDER 


WITH NONLINEAR RESISTANCE 
By 

William Edward Zorumski 
ABSTRACT 

Equations for high intensity acoustic waves are derived from the 
general equations governing a compressible isotropic Newtonian fluid. 

The integral conservation laws of continuum mechanics are used to 
determine a general set of laws, similar to shock wave relations, which 
describe the interaction of these waves with thin porous sheets of 
material. A special case of these laws, where the velocity normal to the 
sheet is continuous and the pressure drop across the sheet is a nonlinear 
function of the normal velocity, is used to study acoustics problems with 
nonlinear material effects. Numerical, approximate, and exact solutions 
are obtained from the one -dimensional problem where two regions are 
coupled through a porous sheet with nonlinear resistance. 

Solutions for scattering from a thin rigid porous elliptic cylin- 
drical shell are found in terms of Mathieu functions. When the acoustic 
resistance of the shell is linear and constsint, the scattering problem 
reduces to an infinite set of linear equations for the coefficients of 
the velocity expansion on the cylinder. Coupling terms in these equa- 
tions are given by integrals involving Mathieu functions which are 
evaluated by exact methods. A special case of variable linear resistance 
is found where the infinite systems uncouple so that closed form 
solutions may be found for the velocity coefficients. Computations are 



made in the intermediate frequency range, where the wavelength is of the 
same order of magnitude as the major axis of the cylinder. Tabulations 
of the scattered and dissipated energies are given for various cylinder 
eccentricities, resistances, incidence angles, and frequencies. A 
perturbation solution is given for scattering from a cylinder with 
slightly nonlinear resistance. This solution is used to obtain the 
scattering of intense sawtooth waves from a cylinder with nonlinear 


resistance. 
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I . INTRODaCTION 


The idealized problem of scattering from a porous cylinder origi- 
nates from the contemporary need for methods to reduce noise from 
aircraft turbofan engines. In reference 1, Marsh has shown that a 
practical method for reducing discrete-tone noise from turbofan 
engines is to Install "broad-band resonators" inside the engine 
nacelle. Because of considerations of weight, safety, and endurance, 
these resonators are usually made of thin porous sheets of material 
(either metallic or fiberglass-plastic) which are fastened to a 
compartmented wall. The cavities behind the porous sheet are usually 
about one-quarter wavelength in depth, since this gives good absorbing 
qualities . 

In general, the greater the exposed area of porous material, 
the more sound is absorbed, so that engine designers must look for 
ways to alter the engine geometry to increase this area. Of course, 
this must be acccxnplished without upsetting the basic flow field 
within the engine, which presumably has already been optimized on a 
performance basis. One approach to this problem has been to install 
thin porous double -walled radial spokes or circumferential rings in 
the engine. The optimum size (from the acoustical viewpoint) of 
these devices is not known, but it is reasonable to assume that their 
dimensions are of the same order of magnitude as the wavelength of the 
dominant tone. In cross section, such a spoke or ring would appear 
roughly as sji ellipse. This suggests that their qualities may be 
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evaluated by studying the tvo-dlmensional problem of scattering from 
a porous elliptical cylinder. 

Measured data In a typlced engine Indicate that the acoustic 
Intensity Is of the order of l6o dB. This corresponds to a pressure 
fluctuation of about one-hundredth of an atmosphere. There Is some 
debate among investigators as to whether this intensity is in 
the llneeu* range (where the classlceil wave equation is valid) or In 
the nonlinear range (where some more exact set of equations must 
be used) . The fact Is that this question may not be answered In 
terms of a statement about acoustic Intensity. Acoustic equations 
are perturbation equations of the general Newtonian fluid equations, 
and each problem must be Judged on Its own merits as to whether It Is 
a linear or nonlinear problem. As evidence of this, severeil Important 
papers on nonlinear acoustics are mentioned here, although no attempt 
Is made to survey this eurea of research. 

In one of the earlier papers (omitting Lord Rayleigh) Flay (ref. 5) 
obtained a Fourier Series solution to the one-dlmenslonal wave 
equation In an unbounded medium. His equations considered a nonllneeu* 
compressibility law (adiabatic) which causes wave peaks to travel 
faster than wave troughs. A viscosity term In the equations limited 
the shape from continued steepening. He concluded that there Is no 
completely stable wave, but that there Is a "neeurly stable” wave 
whose shape changes gradually. In reference 6, Fublnl gives a 
solution to a nonlinear one-dlmenslonal problem which, like Fay's, 
has nonlinear compressibility, but which does not consider viscosity. 
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His solution, which Is given as a series of Bessel functions, 
predicts that a discontinuity will form In the wave. More recently 
Blackstock (ref. k) has generaillzed the work of Fay and Fublnl. His 
work shows that the Fublnl solution Is valid near the source, and 
the Fay solution Is vedld far from the source. Using a weak- shock 
approach, Blackstock has obtained a general solution which Is also 
valid In the tremsltlon region between near and far field. The 
perturbation approach has been used by Coppens and Saunders In 
z*eference 7 to obtain a finite-amplitude standing wave solution In 
a rlgld-walled tube of finite length. Experiments made by them 
confirm theoretical computations for second and third hemDX>nlcs of 
the fundamental wave. Reference 8 Is another example of the use of 
perturbation equations for notillnear equations. In this work, Maslen 
and Moore Investigate strong transverse waves In a circular cylinder. 
The surprising thing about their Investigation Is their prediction 
that strong transverse waves will b? shock free. In contrast to the 
plane wave case where shocks foni. 

Isakovlch, In reference 9> presents a discussion of a 
Sturm-LlovlUe problem obtained from the one- dimensional wave equation. 
He also concludes that second-order terms In these problems will 
always be bounded. 

For the reeuler Intersted In surveying the area of nonlinear 
acoustics, the book by Morse and Ingard (ref. ^), and Beyer's survey 
In Chapter 10 of Mason's "Physical Acoustics" (ref. 2), are 
recosQ&snded as a place to start. 


4 


The nonlinearity which arises from the boundary conditions of the 
problem (rather thun from the governing differential equations) will be 
the focus of attention in the present work. This nonlinearity is a 
characteristic of the porous material. It must be mentioned, however, 
that no time will be spent in detailed consideration of flow and wave 
motion in the porous media. This is a complex subject in itself, but it 
has been developed to a fairly refined state, primarily through the 
efforts of Biot (refs. 10 and ll). In the following discussion, only 
very thin sheets of the porous materials will be considered, so that 
their properties may be characterized by fxmctions which represent the 
changes in flow parameters from one side of the sheet to the other. The 
situation here is analogous to shock wave equations, where relations are 
obtained for computing the discontinuities in flow parameters at the 
shock wave. 

The problem of interaction oi flow with screens or grids is 
important in a number of applications. For example, references 12 and 
15 consider flows •chrough grids of heat exchangers. Screens are used 
for smoothing flows in wind tunnels (refs. l4 and 15 and for 
protecting aircraft engines (ref. l6). These studies were primarily 
concerned with evaluation of a pressure loss coefficient. Reference 17 
is an analytical and experimental study of a strong wave interaction with 
a wire grid. In this research. Center studied a centered rarefaction 
wave as it impinged on the screen. He assumed: 
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(a) quasi- steady flov near the grid 

(b) no heat transfer in the fluid or between the grid and fluid 

(c) effect of grid represented by a drag coefficient which 

depends only on Mach number 

With these assumptions, good agreement was found between theory and 
experiment . 

The theory of dimensional analysis suggests that the drag 
coefficient used by Center shoiild depend on both the Reynolds and the 
Mach number. Careful steady flov investigations by Yates (ref. l8) 
and by Pinker and Herbert (ref. 19) have experimentally conflnned 
this fact. These authors measured coefficients for flows (inside 
the porous screens) varying from laminar, throu^ tiorbulent, to 
choked flow. Both investigations show a dependence on Reynolds’ 
number, but the strongest dependence could be associated with Mach's 
number, which is apparently why Center's assumption was adequate in 
his problem. 

All of the research on screens mentioned above was concerned 
with Mach numbers which are fairly large compared to those found in 
acoustical work where a "large" Mach number would be of the order of 
0.1. Work in acoustics has concentrated on perforated plates, or 
more specifically, on plates with a single perforation, that is, an 
orifice. Also, almost all work in this area has been carried cut 
using iinpedance ideas. These ideas originated in discrete electrical 
ani mechanical systems and are limited to lineeir, steady- state 
problems. In spite of this, the tradition in acoustics has oesn to 
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account for nonlinearity by a modification of the impedance concept, 
which results in a quasl-llnear theory. 

In 1955^ Sivlan conducted an experimental Investigation of 
the Impedance of an orifice (ref. 20). He concluded that the reactance 
of the orifice was only slightly dependent on amplitude (of an 
incident sound wave), but that the resistance Increased with amplitude. 
Sivlan' s work was followed by Ingard in reference 21, who concluded 
that resistance Increases and reactance decreases with amplitude. 

Bles and Wilson made an unsuccessful attempt In 1957 (ref. 22) to 
resolve differences between the work of Sivlan and Ingard. Also in 
1957> Thurston, Hargrove, and Cook (ref. 25) made the first Investi- 
gation of coupling between waves of different frequencies due to 
the nonlinear behavior of an orifice. Their study was of steady 
flow (zero frequency) coupled with a single harmonic wave. 

Reference 2k Is a recent Investigation of the same coupling 
phenomenon by Ingard and I sing. No investigations have been made 
to date of coupling between one or more nonzei*o frequency waves 
striking an orifice or other resistive element with nonlinear 
characteristics. Such investigations would be nearly futile If the 
empirical quasl-llnear Impedance approach were used, because of the 
infinite number of amplitude and frequency combinations which can 
occur In a nonharmonic wave. 

In reference 25, this writer and Parrott presented an acoustical 
theory for thin porous sheets which, like Center's (iref. 17) work, 
used a shock- type relation for the screen. This work was the first 



break with the traditional Impedance approach in acoustics. Experi- 
mental work (ref. 25) confirmed that the theory was useful. In this 
dissertation, a rational derivation of a general acoustic theory for 
thin porous sheets will be given. Solutions for one-dimensional 
problems will be obtained from this. Next, the two-dimensional 
problem of scattering from porous strips and elliptical cylinders 
will be treated. 

The literature on diffraction is extensive. The reader interested 
in this should consult the Important review article by Bouwkamp 
(ref. 26). This survey reviews more than 500 papers published between 
19*4-0 and 195^ • The application of Mathieu functions to diffraction 
problems is discussed in Bouwkamp 's article. Briefly, Mathieu' s 
functions may be used in the study of diffraction (or scattering) 
from elliptic cylinders. His article indicates that research prior 
to 195^ was limited to the studies of the strip or slit (cylinder 
with eccentricity l) . This problem has also been studied with an 
integral equation formulation (ref. 3)* 

Since Bouwkamp 's review, several papers have been published 
which deal with cylinders having eccentricity less than 1. In 1963# 

Yeh discussed the problem of a penetrable strip in terms of a 
Mathieu function series (ref. 27), and Barakat (ref. 28) made a 
study of the elliptic cylinder with various eccentricities. Barakat 
considered the exterior problem using the classical Dlrlchlet 
(sound soft) and Neuman (sound hard) boundary conditions. This work 
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made use of recently computed tables by Baxakat, Houston, and 
Levin (refs. 30 and 31)* These tables are appropriate for 
scattering problems In the lov-to-medlum frequency range. Following 
Barakat's work, Burke and Twersky (refs. 29> 32, and 33) began 
Investigations using low-frequency approximations. 

In reference 29, Burke considered the problem where there Is 
cylinder composed of fluid with a density different frcmi the media 
In which It Is Immersed. Boundary conditions at the surface of the 
elliptic cylinder were continuity of normal velocity and of pressure. 
His paper, like the one by Yeh (ref. 27), considered the coupled 
Interior and exterior cylinder problems. The analysis of this 
problem resulted In an Infinite set of algebraic equations for the 
scattered waves. Burke Indicated the connection between truncating 
these sets of equations and low-frequency approximations. His 
solution gave the far-fleld scattering amplitude to the sixth power 
of the wave number, and the near -field (internal and external) to the 
third power. Terms In the truncated equations could be evaluated 
either from the previously mentioned tables of Barakat, Houston, and 
Levin, from Wiltse and King's tables (refs. 3^ and 35), or frcm 
computations based on known expansions of the Mathleu functions In 
terms of the other tabulated transcendental functions. 

The classical text on the theory of Mathleu functions Is 
McLachlan's work (ref. 36 ). The best Introduction to the subject, 
however. Is probably Arscott's excellent book (ref. 37) • The 
notations are the same In both books, which greatly slopllfles 
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collateral study. Arscott's proofs of the basic theorems and 
McLachlan's profusion of derived equations^ expansions, and 
Identities form a good combination. For the more sophisticated 
reader, Melxner and Schafke’s (ref. 58 ) work Is available, ajid brief 
summaries are given In references 39 f S'Hd 42. Some tables 

are also given In reference 4o, but reference 42, the National Bureau 
of Standards Table, Is the most extensive table available, being an 
extended edition of reference 4l. 



II . DERIVATION OF EQUATIONS AND BOUNDARY CONDITIONS 


2.1 Acoustic Wave Equation 

A derivation of the wave equation may be found In a large number of 
papers and texts (refs. 2 and 3)* The pertinent equations are rederlved 
here In order to clearly define the range of variables for which they 
are valid, and to add continuity to this work. 

Since acoustics Is only a restricted part of the field of fluid 
mechanics, the equations are derived by a straightforward parameter 
perturbation of the general equations governing a conq)resslble 
Isotropic Newtonian fluid. 

The general equations for a klnetlcally euid calorlcally perfect 
gas with constant coefficients of viscosity and conductivity are 
as follows: 

p = pRT (1) 

^+(pqi),i = o (2) 


1 4 . ^ 

15t " 


3 Uji 


(5) 




[p(^ . Cpl)] * [p(^ . Cpl)4i],i = I 


- I p[4k,k1i] ,1 + p[(ll,J + + KT,ii 


(M 


10 
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since only small fluctuations of the fluid state are to be 
considered, the ambient pressure, density, and temperature, 

Pq, Pq, and Tq, respectively, will be used as reference quantities. 

A reference length and time are also needed, and for these the 
reciprocals of clrculeu: frequency o), and wave number, k, will be 
used. 

Thus, denoting dimensional variables by ('^)'s, and dimensionless 
ones by (x)'s, the following set of dimensionless variables Is 
Introduced: 


II 

X 

(5a) 

t* =1 

(5b) 

CD 

^ 5i 

(5c) 

P* = PoP 

(5d) 

P* = pQP 

(5e) 

T* = TqT 

(5f) 


A few remarks about the above quantities are needed here. The 
frequency range of Interest Is of the order of 1(P Hz, perhaps 
100 Hz to 10,000 Hz. Correspondingly, the wanre lengths will vaiy 
from. 55 cm to_0.55 cm and the wave numbers from 0.2 cm“^ to 20 cm"^. 

will be roughly the speed of sound, or about 550 meters/sec. 
Thus, ^ Is a Mach's number, and will be quite small, p Is the 
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pressure fluctuation in atmospheres, which is also assumed to be very 
small . 

The reference pressure, density, and temperature will be chosen 
to satisfy the perfect gas law. 


Po " Pq^o 


( 6 ) 


The dimensional equations will now be converted to dimensionless 
variables. Equation (l), with the condition of equation (6), is simply 


p = p T 


( 7 ) 


The continuity equation is of the same form as before, namely 

9p I ^ 0 

^t 


( 8 ) 


Equation (3)^ the momentum equation, contains two parameters 
after the changes of variable. 

- /k^ Sp , ^ \ 

Dt 1 3 ^xj dxi 

The energy equation also brings in two additional parameters. 
After transformation it is 


( 9 ) 


^t 


-M . 


o )\2 








+ p +-T^|qi 


Po 8p 




Po 


(^) 


2 

5 SX4 


- 






^^xj dx 




»iKT, 


S T 


Po(^)^ ^1 8xi 


( 10 ) 
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Consider the set of parameters which appear in the transformed 
equations (7) through (lO). They are 


Po Ku ^^o 


and 




The first term, 



is approximately 


1 
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This is a first order term. The second is an inverted Reynolds number. 



This number is of the order 10”^ for the wavelengths of interest 
in air, so that it may justifiably be called a small quantity. 


is roughly 3 for air, so it must be considered a first-order 


CpTp 

{if 

term. — ; — - varies from about O.O 5 x 10"^ to 3 x 10”^, the larger 
'h)\3 


<=o(^) 


value corresponding to the larger wave number, or frequency, thus it 
may be considered to be a small quantity. Using the dimensionless 
groups discussed above, the following quantities are defined: 






cK 


(lla,b) 


(11c, d) 



Ik 


The new quantities defined In equations (ll) may now be thought 



of as of order 1. If c Is not chosen smaller than 10 , will 

be roughly 1 or greater, auid k will be roughly less than 3 . 

Asymptotic solutions for equations (7) through (lO) may be 
found In the form 


p ~ 1 + + €^2 + 


- _ 2 — 

P ~ 1 + + 6 p2 + . . 


T ~ 1 + CT, + + . 


qjL - 0 + €q^^ + + 


( 12 a) 

( 12 b) 

( 12 c) 

(I2d) 


Substituting these asyiq>totlc expansions Into the governing 
equations gives the desired sets of first and second order equations. 
This Is a routine procedure except for equation (lO), the energy 
equation. Here, It must be noted that the leading term In the 
ten^erature expansion gives a large reference energy level. Because 
of the continuity equation, this term drops out of the energy 
equation. This can be accounted for by substituting CT^^ + c^T 2 Into 
the energy equation, that Is, by using a relative teoq)erature . 

If this Is not done, the resulting prediction would be that the 
speed of sovind Is / FT, Instead of the correct value of \J 7 HT. 

Keeping this In mind, the govenilng first and second order asyiqptotlc 
equations may be found. They are as follows: 

Pi ~ (Pi + = 0 


(15a) 



The first order set of equations ( 15 ) may now be manipulated to 

obtain wave equations. 

Since C_7 = — — 

P 7-1 

the equation 

^ bt 

is a direct result of the perfect gas law and the energy equF.tion. 
With this, we may find the wave equation 
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2.2 Porous Sheet Equations 

A form for a set of equations which relate pressure, velocity, and 
temperature (or density) on opposite sides of the thin porous sheet is 
suggested by the integral form of the laws of mechaziics. The integral 
mass, momentum, and energy equations may be written in terms of the 
dimensionless variables in equations 2.1(5) as 



Consider a control volume which encloses the fluid flowing inside an 
element of the porous sheet, as shown in figure 1. The dimensionless 
thickness of the sheet, h is much less than 1. The control volume 
is fixed with respect to the sheet, which may be moving with a 
velocity, Since the control volume encloses only the flviid, and 

not the solid, the surface integrals extend over the interior 
wetted surface of the solid, and over the faces of the element ^ere 
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fluid is passing through. The velocity of the fluid may be 
expressed in terms of a relative velocity as 

^i = Vi + v^ (i+) 

Since the flow parameters will vary greatly from point to point 
within the porous element, they will be broken into average (over 
an area) and fluctuating parts. The average of a quantity is 
defined as 

^ = lJqdA ( 5 ) 

where A is an area whose size is small cosq>ared to a characteristic 
length of the problem, but large enough to give a meaningful average. 

The fluctuation of a quantity is then 

Q' = Q - S (6) 

In this problem the Integrals will be taken over the faces of the 
control volume shown in figure 1. 

The volume Integrals in equations (l), (2), and (5) give time rates- 
of -change of mass and momentum and energy stored inside the control 

volume. These integrals and the surface Integrals on the right side 
of the equations cannot be evaluated analytically. About all that 
can be said is that they are of two types: One with a time derivative, 

and one without. 

The surface Integral in equation (1) may be written as 
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f_ Pn^q^ dS = f priiVidS+ h“^|p(Vi+ v^) + (p' V;[)] -[p(Vi + Vq^)+ (p'V^ | 
"S ^Sint h -^oJ 


+ similar terms for faces 2 and 3 of element 


Consequently, the mass equation can be represented by the form 


( 7 ) 


^[p(Vi + v^)J = h/^[p(V2 + V2)J + ^ 


+ Mr 


(8) 


where 


a[p(Vi + = p(Vi + vi)j - p(Vi + V3^)|_ 


(9) 


o 'h 

The first group of terms or the right of equation (8) represents 
flows in the plane of the sheet. The time derivative arises from the 
previously mentioned volirme integral, and the last term, M2, includes 
the integral over the internal wetted surface and terms like (p*v'). 
Similar forms are suggested for the momentum and energy equations, namely 


f~ 

A p(Vi + V2^)(Vj^ + ^ ^ + l2i 

L 7 J ^ 

+ h^^[3(V2 + V2)(Vi + v^) + =1 + — [^(V^ + v,)(Vi + Vj^) + =1 

BX2L 7J 8xjL 7J 


( 10 ) 


and 
A 


+ + vg) ((y . t .l f . + c^t) + (V 2 + V2)|j 

+ 2 ^ * ^''5 + S'ljj 


( 11 ) 
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The functions M, L, and E must be evaluated experimentally 
In terms of the average flow parameters. Once this Is done, 
equations (8), (lO), and (ll) give the necessary relations between 
flow parameters on either side of the porous sheet. Since the 
equations involve the velocity of the porous sheet, this quantity 
must either be specified or accounted for by adding appropriate 
equations for the motion of the sheet, such as plate or shell 
equations . 

For the balance of this study, the porous sheet will be 
considered to be fixed. Cross-flow terms will be neglected since 
the sheet thickness is small. 

With these conditions the sheet equations become 

^[^1 _ ^ + Mp ( 12 ) 

at 


and 



(15) 

(Ik) 


The average flow variables at each face of the porous sheet are 
assumed to equal the acoustic variables, that is 


P = 


(15a) 


(15b) 


$ = T, 


(15c) 


Vi = q^ 


(I5d) 


and 
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Using these equations and the asymptotic expansions 2.l(l2), 
equations (12), (l5)# and (lU) become 


« 4«ll] = ^ + Mj, + 0(e2) 


^t 


8ii€ A 


fpil ^ ^ ^ L + 0(c2) 
ot 


L^J 


g2 m 


‘111"" “ 

_ 7 ot 


(16) 

(IT) 

(18) 


A limited amount of experimental work (ref. 25) indicates that 

p 

Ml and M2 are of order e . The energy terms can only be of order 
and, since the left side of (17) has the factor L |-[ and L21 

must be of order e, while Li2> Li3> ^22^ order e^. 

I'll found to be proportional to q -| 1 , and L21 is a fianction of 

Mach number with a slight dependence on Reynolds number ^ich can be 
relegated to the second order spot. Consequently, for acoustics 
problems involving thin porous sheets of material, the bo\indary 
conditions at the material become 


Apii] = 0 


(19) 


€ 



- 

em — ii + AP 
8t 


pii] 


(20) 


Here, m is a constant and AP is a nonlinear function. 



III. ONE-DIMENSIONAL PROBLEM 
3*1 Pressure Drop Function 

Little has been said to this point about the origin of the 
parameter e. It could be associated with a forcing function, but in 
riany acoustics problems, these do not appear explicitly. In this paper, 
it will be associated with a nominal sound pressure level. 

The sound pressure level, or intensity, in decibels, is defined as 



( 1 ) 


( 2 ) 


SPL = 20 Log]_Q(lO^° ep;^) " 5-0 


In order to make the dimensionless pressure vary from 0 to 1, € 
may be chosen as some convenient power of 10 to give the correct SPL 
range. Table I shows the correspondence between SPL and e when the 
dimensionless pressure is one. 

The nonlinear pressure drop function in equation 2.2 (20) depends 


on €. If that equation is divided by €, all terms must be of order 1. 


The pressure drop term is 



Figure 2 shows a typical plot of 


this function for several values of €. Each curve shown may be 
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approximated fairly well by a straight line. This does not mean that 
the nonlinearity is small. The slopes of the approximating lines in 
figure 2 vary from about 0.2 to J.O so, if this slope is thought of as 
an effective acoustic resistance, it must be concluded that the variation 
in resistance is large. Nevertheless, the fact that some straight line 
approximation is possible may permit perturbation methods to be used. 


5.2 General Equations for Two Coupled Regions 
In equations 2.1 ((17) and (l8)) and 2.2 ((19) and (20)), let 


and 





m 

— 

7 


(la) 


(ic) 



(Id) 


The one -dimensional wave equation is then satisfied by 


U= ^ (2) 

n=-<» 

oo 

Z , -Inx inxv Int 

(u^e - v^e )e (5) 

n=^ 

In equations (2) and (5)# and v^ are complex constants with 

(Ua) 


» (u.n) 
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V 


n “ 



(4b) 


Each term with coefficient represents a right -moving harmonic 
wave, while the terms give left-moving waves. 

Consider the two regions shown in figure In each region there 
may be a periodic solution which is given by equations (2) and (3). At 
the sheet, there are the conditions 

Ul = Ug (5) 


_ _ . , AiCuJ 

Pi - - Pa 5^ + 


( 6 ) 


and at the extreme ends of the regions there are general impedance -type 
boundary conditions 


Pin " ^In^ln ''' ^In 

(7a) 

P2n “ ^2n^2n * ^2n 

(7b) 


Here the subscripts 1 and 2 indicate the region while the subscript 
n indicates the harmonic. 

If equations (2) and (3) are substituted into (5) through (7) and 
use is made of the orthogonsLLity of con 5 )lex exponentials, there resiLl.ts 

(8a) 


“In * ’'in “ “2n + ’'2n 


(8b) 




Figure 3.- One-diiLensional problem for two regions coupled throiogh 

a nonlinear porous sheet. 
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- ’'In) - (>^2n ’ ''2a'l 



-int 


e 







.imt 


dt 


+ inPA(uj^ + Vj^) 


(8c) 


-an--'”" - W'”" = 22nCw-^® - V2n«^“®] " P2n 


(8d) 


5.5 Approximate Methods of Solution 
T3ie above equations are an Infinite set of nonlinear algebraic 
equations. They may be specialized to simulate a variety of situations. 
For example, If = 0, and Z 2 J. = 1# the region 2 will contain only 

right -moving waves. This simulates an Infinite region where no waves 
are reflected. If Z 2 n Is Infinite at X 2 = it/2, a rigid wall Is 
simulated. The solution to these equations may be approximated In a 
variety of ways. 

They are difficult to work with numerically, however, because In 
any Iteration process, the Integral term must be evaluated a large 
number of times as the Iteration proceeds. In spite of this, the 
following technique has been successfully used to obtain numerical 
solutions. 


3 . 5-1 Steepest Descent Iteration (ref. 43) 

If an Initial estimate for the solution Is made, that Is, 

’’■■'n# ^2n# specified, then equations 3*2 ( 0 ) may be used 

to compute the nonhomogeneous terms, P-i^,, Pp^. The problem solution 

Is then given by the minimum of the function 

00 

n=-» 


( 1 ) 
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A steepest descent technique was used for locating the mininiuin of 
equation (2C). This method amounts to following the gradient of the 
potential function, S, during the iteration. The reader interested in 
further details of the method may find them in reference 


3«5*2 Collocation ^fethod 

The integral which causes numerical difficulty may be eliminated 
by adopting a different criterion (ref. for satisfying the pressure 
drop condition. 

For brevity, we introduce the notations 


= ^In + ^In = ^2n ^2n (l) 

and 



In terms of these variables, eqxiatlons 3*2 (8) become 


^Pn = Vn + Pn 


( 5 ) 


wnere 


Z 

n 


Zin cos(nx^) + i sin(nX]L) 
cos(nx2^) + iZ]_j^ sin(nxj_) 


^2n cos(nx2) + i sin(nx2) 
cos(nx2) + IZpy, sin(nx2) 


- 


(M 


and 
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^In ^2n , , 

p (5) 

^ cos(nX]^) + iZQjj sin(iiX]^) cos(nx2) + iZ2n sin(nx2) 

In equatJ.on (5)# the are the complex Fourier coefficients of 

the steady flow pressure drop, Note that the inertia term from 

equation 2.2 ( 20 ) is included on the right of equation ( 5 )* To satisfy 
the pressure drop law continuously would require that 


Mt) 


fiP(€V(t)} 

e 


( 6 ) 


An approximation to the condition (6) is to use a least-squares 
criterion at discrete points in time, that is, use the minimum of the 
function 



( 7 ) 


to approximate the problem solution. 

The instantaneous pressure and velocity are given in terms of their 
Fourier coefficients by 

L L 

= ^ ( 8 ), ( 9 ) 

k=-L k=-L 



where 


0 < L < 2 
2 
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Evenly spaced intervals are taken in (7), so that 

_ 2jtn 


tr, = 


N 


so that the potential function, S becomes 


N^l L 2^ 

n=0 U=-L 



) ( 11 ) 


A powerful computational algorithm, known as the Fast Fourier 
Transform, has been developed for working with the sums which appear 
inside equation (ll). Details of this are available in references 45 
through 48. 

Simeonov's method (ref. 53) or the steepest descent procedure could 
be used to minimize the function given in equation (ll) with respect to 
the Fourier coefficients, Vq. Vflien this is done, however, a summation, 
which is analogous to the integral in equation 3*2 (8c), must be made 
repeatedly during the iteration process. This may be avoided by using 
equation (6) to solve for the velocity at N discrete i>oints in time, 

■^n* 

The inversion of equation (9) gives 


N-1 

Vn = V(tje‘ 
m=0 


2jrinm 


N 


( 12 ) 
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Using this result with equations (3), (6), and (8) gives 


N-1 


^ Ap(eV(t )} 

* Htn) 


( 15 ) 


m=i 


where 


L 2 


■'nm 


21 


jtik( n-m) 




(1*1) 


k=-L 


and 


L, 2jtikn 
N 


P(tn) = 2. 

k=-L 


(15) 


The solution of the non l inear algebraic equations (13) satisfies (6) 
at a finite set of times. To satisfy the pressure drop relation con- 
tinuously would require an infinite set of equations. The advantage of 
( 13 ) is that the transform of the nonlinear function does not appear. 

A variety of techniques may be used to solve (13)* Since it gives the 
minimum of (7), which is zero, the steepest descent procedure, or 
Simeonov*s method of descent could be used, or a technique of successive 
approximations might be a feasible method of solution. A disadvantage 
of this set of equations is the need to invert the array G^, which 
may be a very large order matrix. 


3.3‘3 Linearized Perturbation Equations 
The data shown in figure 2 indicate that the pressure drop function 


may be written as 


55 


— i — = + sifq J 


R, r, and & each depend on the parameter e. For purposes of 
analysis, 6 is taken to he a small parameter which is large when 
compared to £ hut smell compared to one. 


6 = < < 1 


R = 0(1) 


r&iJ = 0(1) 


Using equations 3*2 (8) again, and assuming a power series solution 


in terms of the parameter 5, 


( 0 ). . ( 1 ) 


^In ~ ^In + S^ln + • • • 


V ~ . . . 

In In In 


"2n ~ • • • 




gives the zeroth and first-order linear power series equations 
(O) -inxi ^inxi _ ["(O) -inxi Jo) inxi"! 

^In® ^In' ^InPln® ^ ^In® J ^ 

(0) (o: (0) 

^In ^ ^In - ^2n ^ ^2n 


(4h) 
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(4? - ''Is’) - (4S’ - 4S’) = (R + i-“0a)(4S’ + 4S’) ('“=) 

40 )^-lnx 2 . 40 ),lnX 2 ^ 




+ 4i’ = 4i’ + 4S’ 


<4n^ - 4n4 - - ’'aj4 = (R + lnp^)(uW + v(l)) 


2lt./.* 


g-intri 


•[v(°)(t)] 

,U),-lnx2 . , Z^[u^^)e-i-2 . v^^)ei-2] 


dt 


(5a) 

(5b) 


(5c) 


(5d) 


The advantage in this approach is apparent. Solving each set of 

linear equations is a trivial matter. The integral term which causes 

computational difficulty in an iteration scheme must be evaluated only 

once to obtain the nonhomogeneous part of the first-order equations. 

There is a danger in using this method in that the choice of € will not 

give a solution such that qj^2_ varies as expected. For example, in 

figure 2, when e = 10” it is assimed that q-j^j. will have a maximum 

value of 1.0. When e = 10 , (Q-tt) should be about 0.4. It is 

XX tnsjx 

possible to extend this reasoning to obtain an iteration scheme based 
on power series equations. 

3.4 Exact Solutions 

In terms of the notation introduced in equation 3.2 (l), the 
one -dimensional wave equation is 
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^ ^ 


and 




( 1 ) 


( 2 ) 


The classical solution (ref. 49) to (43) is 


U = f(t - x) + g(t + x) 


( 3 ) 


Figure 4a depicts a problem where an incident wave, f]_, is reflected 
from a porous sheet which is backed up by a perfect absorber. At the 
porous sheet, we use a special case of the pressure drop law, with 
Pa = 0. 


Pi " 





The acoustic velocity is continuous at the sample and there is only 
a right moving wave behind the sample. 


fi(t) + gj^(t) - fgCt) = 0 


fl(t) - gi(t) - f2(t) = 


AP[e(fj^(t) + gj^(t))j 


(5) 


( 6 ) 


Let 


p (u) = + 21J, 


€ 


(7) 
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then, by subtracting (5) from (6) and adding 2 f 2 ^(t) 
the resulting equation, it may be shown that 

to each side of 

Ps(fl(t) + Bi(t)) =2fj^(t) 

(8) 

8l(t) = - f^(t) + 3"^(2fj^(t)) 

(9) 


The exact solution, given by equation (9), depends on the possibility 

of inverting the function PgCu). This presents no practical difficulty. 
^1 £ttJ 

Since — ^ — * must be constructed from experimental data, there can be no 
loss in accuracy by directly constructing the Inverse function, 3^^, 
from the same data. 

For the problem depicted in figure 4b, the equations which govern 
the transmitted and reflected waves are 


fl(t) + gi(t) - fg(t) - ggCt) = 0 


fg(t - w) + gg(t + w) = 0 


fl(t) - gi(t) - fgCt) + ggCt) = 


APf (fl(t) + 6i(t))] 


( 10 ) 

( 11 ) 


( 12 ) 


The incident wave, f, is considered to be given. Equations (lO), (ll), 
Euad (12) may be used to eliminate f 2 and g 2 to give an equation 
for g^. First adding and then subtracting equations (lo) and (12) gives 


fg(t) = fj^(t) 


1 £j{|e(fi(t) + gi(t))] 

2 e 




^ ^€(f^(t) + g^(t))] 
2 


(15) 

(iM 


e 
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Substituting these equations Into (ll) gives 


^ AP0(f3^(t - w) + gj^(t - w))] 


f, (t - w) + g, (t + w) - - 

2 

^ Z^0e(f^(t + w) + + w))|] 


+ = 

2 


0 (15) 


Since (15) holds for all values of t, it must hold at t - w and 
t + w, thus 


fl(t) + gi(t + 2w) - I 


^ ^[e(fl(t) + gi(t))] 


T ^[e(fi(t + 2w) + gi(t + 2W))"] 

+ ± : 1 d = 0 ( 16 ) 

2 e 


, Z^<i(fl(t - 2w) + g^(t - 2w))1 

f^(t - 2w) + g^(t) - i i i 

2 e 

1 Z^[€(fi(t) + gi(t))1 

+ = = 0 (17) 


Adding (l6) and (17) gives the identity 

^l(t) + g]^(t) + fj^(t - 2w) + g^(t + 2w) = 0 ( 18 ) 

Let f^(t) have a period equal to 4w, then 

^l(t) + gi(t) = - (f^(t + 2w) + g^(t + 2w)) (19) 

Now, if is odd, then substituting (19) into (l6) gives 
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fl(t) + gi(t + 2w) = 


AP[£(fi(t) + gi(t))] 


( 20 ) 


Letting t go to t - 2w in (19) and substituting into (17) gives 


f^(t + 2w) + gj^(t) = - 




(21) 


Adding and subtracting f]^(t) in (21) gives 


/ff[6(fi(t) + gl(t))] 


+ (fj^(t) + gj^(t)) = fj^(t) - fj^(t + 2w) (22) 


Now, if 


Pl(u) 




(23) 


the exact solution for the reflected wave Is 

gl(t) . - f^(t) + Pi^(fi(t) - fi(t + 2w)] (24) 

In the geneiel case, where the backing is not a perfect absorber, 
or where the rigid wall is not exactly l/4 wavelengths away, the waves 
behind the san5)le are related by the convolution Integral 

n 

0 (t - T)fg(T)dT . G*fg (25) 

n 



where 


G(t - T) . i 


y / 1 - ^(0)\ in(t-T) 

n4l-oU + 


(26) 
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and Zjj(o) is the acoustic impedance in the backing cavity, evaluated 
at the sample. Thus, in general, the problems of wave reflection from 
a porous absorber with nonlinear resistance depends on the solution 
of the set of nonlinear functional equations (lO), (12), and (25). No 
general theory exists for the solution of functional equations, although 
a number of examples are discussed in the book by Aczel (ref. 50). 

A numerical approach to their solution gives equations like 3»3.2 (13) # 
that is, large sets of nonlinear algebraic equations. 



IV. SCATTERING FROM A POROUS ELLIPTICAL CYLINDER 


4.1. Coordinates, Governing Eqmtlons, and Boundary Conditions 
Flgvire 5 shows the elliptical cylinder coordinates which will be 
used in this section. These coordinates are define by 

X = h cosh 5 cos n (l) 

y = h sinh | sin t) (2) 


As shown in flRure 5^ curves with i = constant are ellipses with focal 
points at x = ± h. Curves \^ere t] = constant give hyperbolas with 
focal points which correspond to those of the ellipses. 

Let = n^q.^ and denote a unit normal to the elliptical cylinder 
by n^. Then the governing equations are 


^^p _ 


(5) 


^ 

A[niqij = 0 

An ^(n^q.) Ap[en-qJ 

= "A + -T^ 


(M 

(:?) 

( 6 ) 


When these conditions are expressed in elliptical coordinates, 
they become 

= 2k^(cosh 2| - cos (7) 

3t, 2 at® 

^ ^ ( 8 ) 

^ \/2k^(cosh 25 - cos 




Flgxire 5»- Elliptical coordinates. 




Here, is the velocity along a curve q • constant, denotes 
the surface of the porous cylinder, and 


k = h 
2 


k,2. Derivation of Solutions to Mathleu's Equation 
We seek a periodic solution to equation 4,1(7) of the form 


00 

Pn(^^T) 


For these solutions the wave equation becomes 


— ^ ^ + (nh)^(cosh2 i - cos^ ri)p (| t]) = 0 

and the velocity Is related to the pressure by 


1 2k (cosh 2| - cos 2y\) ri=-oo 


in 


y ^ 

Z_, in 


T I T ' U in 

\/ 2k (cosh 2i - cos 2r\) n =-oo 

The slngulau: terms In equations ( 3 ) and ( 4 ) require that 
Pq(I»^) “ constant. Although these terns add a small circulation 
and divergence to the flow field, they will be omitted for the 


( 12 ) 

( 1 ) 

( 2 ) 

(5) 

(4) 


remainder of this discussion. 
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Sepeui^ble solutions of equation (2) are of the form 

Pn(^>n) = (5) 

Inserting this form Into the wave equation gives 



0n*vT) + (t - Sjj cos^ n) = 0 

(6) 

and 

- (b - Sjj cosh^ 0 ’^n(^) = ° 

(7) 

where 

o r,2v2 

— n h 

(8) 


For the development of the theory of Kathleu functions, the 
notation of the National Bureau of Standards Tables, reference h2, 
vlU be used since later computations are to be based on these tables. 
The method of derivation Is t hat which Arscott gives In his book 
(ref. 37). 


k,2,l Circumferential Fourier Series Solutions 
In equations ^.2(6) and 4.2(7)> which are called Mathleus 
equation and Mathleu's modified equation, b Is a separation 
constant which depends on s^. For to be periodic, b will 

take a set of values. To each of these values, there corresponds 
an Independent solution. These solutions fom a complete orthogonal 
set. 

Arscott has given a clear and rigorous proof that these solutions 
are of the fora 

00 

Se2r(sn,n) cos 2kri ; 

k=0 


( 1 ) 


1;5 

=£ °4k+l^^ cos(2k+l)Ti ; 
k=0 

00 

So2r+i(sn.Ti) sin(2k+l)Ti ; b02j,+i(sn) (5) 

k=0 
00 

S 02 j.+ 2 (sn^n) ^°2k+2^^ sln(2k+2)ii ; b 02 j.+ 2 (sn) (^) 

k=0 

In equation 4.2.1(14-) ve have deviated slightly from the NBS 
notation by using the subscript (2r + 2) Instead of 2r. This 
allows the suimnation to extend from o to « In all cases. 

Functions and coefficients with the same numerical value are still 
as defined In the MBS tables. Note that the even functions Include 
the letter ”e” In their symbols^ and the odd functions use the 
letter "o.” Methods of evaluating the coefficients De, Do, and 
separation constants, be, bo are given in the KBS tables. The 
coefficients axe normalized according to the condition 


Sej,(sn,o) = 1 

dSo, 


dT] 


= 1 


T1=0 


This normalisation gives the orthogonality conditions 

p2jt 

/ Sei.(sn,T]) SOm(sn,Ti)dii = 0 


o 




Ser(sj^,n) Sem(sn,Tl)dTi 



SOr(sn,Ti) SOi^iC n)dTi = 



m / r 
m = r 
m ^ r 
m = r 


(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 
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4.2.2 Radial Bessel Function Product Solutions 
To generate the radial solutions for the modified Mathleu 
equation, use may be made of the formula 

r2i( 

♦n(5) = Cn / Gn(6.l) 0n(’l)d'l (D 

^ o 

It may be shown that the function given by equation (l) satisfies the 
modified Mathleu equation If Is any solution of the wave 

equation and 0 q ( tj ) Is any of the periodic solutions of Mathleu *s 
equation . 

Appropriate kernels for equation (l) are the cylindrical waves. 

To get these, the wave equation Is transformed Into i>olar coordinates 
where It Is 


1 i. (r + 

^ 3r ' dr / 


d^p(r,9) 


+ n2Gn(r,e) 


= 0 


( 2 ) 


where 


X = r cos 0 
y = r sin 0 

Equation (31) has the solutions 


(5) 

(M 


/ -1 o \ COS 

Gn(r,0) = 7^0 ; 7^ = 0,1,2,... 


(5) 


In reference 51 ^ Watson gives a proof of Graf's ewidltlon theorem 
for Bessel functions, which Is ^ 


Ha(“) Sin {■“«'> 


( 6 ) 



In this formula 



03 = 

- 2Zz COB 0 

(7) 


03 

cos t = Z - z cos 0 

(8) 



03 sin ^ = z sin y 

(9) 

If we let 


0 = 2T) 

(10) 



z = S gt 
2 

(11) 



z = - 

2 

(12) 

then 


03 = nr 

(13) 

and 

COS 0 = 

cos T] cos - sin r, sin i|f 

(14) 


sin 0 = 

sin T sin ^ + cos ^ cos \|f 

(15) 


sln(^'*i} =£ (-l)X-Hn(^ e"^) sin 


Equations (6) and (l6) hold for Hankel functions of the first or the 
second kind. 

Although euiy solution of the wave equation serves as a kernel 
for equation (l), it is possible to obtain trivial solutitns. 



i+8 


Nontrivial solutions are found by selecting kernels which have the 
same period as the generating solution, and which are even or odd 
as is even or odd. 

Corresponding to each of the circumferential solutions, ^.2.l(l) 
through ^.2,l(l4-), there are the radial solutions 
p2jt 

He2r(sn,5) = C2r / HQ(nr)Se2j.(s^, Tj)dTi; be2j.(sn) (l?) 

^o 

p2jt 

He2r+l(sn^O ^2r+l / cos 0 Se2j.^.]_(sj|^,l)dilj be2p^3_(sn) (l8) 

o 

p2n 

H 02 r+i(sn,l) = ^ 2^+1 / Hi(nr) sin 0 So 2 r+i(sn,^)<n; l302r+l(sn) (19) 

^o 

2jt 

H°2r+2(®n>0 = d2r+2 / H2(nr) sin 20 So2j.+2 (sn^ b 02j.+2(sn) (20) 

o 


In these equations, Hankel fimctions of the first or second kinds 
may be used, so that the complex conjugates of these expressions are 
independent radial solutions. By using the trigonometric identities (14) 
and (15), Graf's addition theorem (16), and the expansions for the 
periodic Mathieu functions, the radial solutions are reduced to the series 


He,,(s„,0 =1^^ I e-A; .e,r(s„) 

r\ ^ ' 


° m=0 


( 21 ) 


1 m=0 L 


+ 




; be 


2r+l 



( 22 ) 
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Ho 


2 r+l 


(sn;0 = 



(-iTdo 


(2r+l) 

2 m+l 







^m+l 




j ^°2r+l^®n) 




(-1)^ 

Doo 


m=0 


(2r+2) 

2 m+2 






^°2r+2^®n^ 


( 25 ) 


(24) 


Equations (21) through (24) are the railal solution which will be 
used for later computations. Since the comple-c conjugates of those 
solutions are also solutions, these equations each give two 
Independent radial solutions. To distinguish these we add super- 
scripts such as He^^] ^®2r* 

Asymptotic expansions for the radial solutions with large 
arguments may be found from the known expansions for Hankel functions. 

In reference 52, Hildebrand gives the formula 

± i fx . 

h(1^2)(jj) e L ^ J (25) 

P Vnx 

With this formula and the values of the Bessel function at 
the origin 


Jp(o) 



p = 0 
V ^ 0 


the asymptotic formulae for radial solutions may be shown to be 


( 26 ) 
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. A^‘ t 

rt 

(27) 

(1,2), f 1 \r t '■P? 

eS 

5 «" 
■ U _ 

( 28 ) 

(1,2) (-!)'■ 

^°2r+l <-== — - ® ^ 

\|^ 


(29) 

(12) ( i^r+1 ± i 1 ^^ 

Ho,J ’ '(s,,|)=l^e l2 

-e] 

■'( 50 ) 

Derivatives of the radial solutions may be found by differentiating 

equations (2l) through (2^4-) and using the Bessel Function identities 

for derivatives 



H^(x) = Hp.i(x) - E Hp(x) 


(31) 

H^(x) = -Hp+i(x) + E Hp(x) 


(32) 


The general separable solution for may no\^' be given as 

00 

P„(6,l) Ser(s„,il) 

+ |b4>Ho^^’(s„,{) + SOr(6n,T|)j (55) 

Pn(5,l) and P(_n)(5,^) must be conjugate fimctions in order for 
r(5,1,t) to be real. This condition introduces the following constraints 
on the undetermined coefficii.nts of equation (55) 
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a(i) 

-nr 


' nr 



B 


( 1 ) 

-nr 





' nr ’ 


' nr 


) n = 1,2,3 


J 


(3^a) 


r5Ub) 


(3^c) 


i^hd) 


Equations (5^) show that it is sufficient to determine the complex 
coefficients for positive values of n only. 

Substituting (53) into ^.2(l) and using the asymptotic fomiula for 
the radial Mathieu fvinctions gives a far-field pressure formula 


pU,n,t) 



i|n(ke^+t)-0j 

e 


+ A 


( 2 ) 


nr 


+ B 


( 2 ) 


nr 



i^(ke^+t )-9j.J 
e 


(35) 


wnere 


9 


r 



r even 
r odd 


This equation clearly shows that the coefficients and B„^ must 

^ nr nr 

be associated with diverging waves, while A^^^ and are 

associated with converging waves. 


52 


It is often convenient to work with the real and imaginary parts of 
the radial solutionc. In reference 42, these are defined as 


ne(^^(Sn,0 = Jej.(sjj,0 + iNej.(sn,l) 


( 56 ) 


Ho^^^Sn,!) = JOj.(sn,0 + iNOr(sn,0 (37) 


The Wronskian of these solutions may he shown to be equal to 1. 
Therefore 



JejJJe^ 

- Ne^Je^ = 1 

( 38 ) 

and 





JOj-NOr 

- NOpJor = 1 

( 39 ) 

for all values of 

Sn> 

r, aind 

(40) 


Finally, it may be seen from equations (2l) to (2k) that Je^. 
and JOj. are even and odd functions, respectively. Ne^. and NOj., 
however, are neither even nor odd. 


4.3 Interior Cylinder Problem 
The pressxire within the cylinder is given by substituting 
4.2.2(55) into 4.2(1) 

oc oc 

n=-oo r=0^ 

+ + B^)Ho(2r|sorJ (l) 


0 



55 


In order for the pressure and its gradient to be continuous on 
the x-aocis (5 = O), 


p^(0,T),t) = p^(0,-Ti,t) 

and 


( 2 ) 

( 3 ) 


By using the orthogonality equations of the complex exponential and of 
the circumferential Mathieu functions, conditions (2) and ( 5 ) may be 
shown to be equivalent to 


4 ^^ 

= a(2) 

nr 

(4) 

nr 

nr 

(5) 


Let the velocity normal to the surface of the cylinder be given 
by the expansion 


00 00 


I I Ser(T]) + VnrSOr(T]))« 


int 




_ n=-<” r 


=0 


( 6 ) 




2k^(cosh 2lo - cos 2r\) 


If equations 4.2(3) and (l) are used to compute the velocity at 
the surface and the result is equated to the expression in (6), then 
two relations are obtained for the and coefficients. 
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(7) 


2B^4o^'(sn,lo) = -i“V 


(8) 

By using equations (4), (5), (7), and (8) to eliminate the 

coefficients 

Anr ®nr 

from equation (l), we obtain 




00 00 



= - 

y InV |c4.(|)Unj.Ser(>l) + 

int 

e 

(9) 


n=-i» r=0 J 



where 




and 

i Jej.(£) 


(10) 


q 1 _ Jor(l) 
Pnr f . . 

So) 


(11) 


4.4 Exterior Cylinder Problem 



Consider a 

periodic wave 'vdiich is moving in the positive sense 


along the line 

X + iy = R e^®* 


(1) 


This wave may be expressed as 

00 

p(5,1,t) = ^ Pn (2) 

n=-oc 


The normal distance from the origin to a point on the wave Is given by 
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R = X cos a + y sin a 


( 3 ) 


or 


nR = y/s^(cosh I cos t] cos a + sinh 5 sin T sin a) (^4-) 

It is desired to represent this plane wave in the form of 
equation (65), thus we write the tentative expansion 


-inr 


r=0 


Sej. + + Fo^^^Ho^^^ So^ (5) 

nrr nrr ^ nrr r r rf'-" 


The orthogonality conditions on the periodic Mathleu functions 
yield simultaneous equations 


Nr [j 




'] -r 


-inR 


e Sej,(T))dT] 


( 6 ) 


“r e'^'^SOr(ri)cn 


(7) 


The exponential term in equations (79) ®nd (80) satisfies the 
wave equation, and since R is invariant when a and tj are Interchanged 


.2 „-inR ^-inR 


Therefore, the integral terra in (6) satisfies both Mathieu equations 


( 8 ) 


2 

+ (bej, - Sj^ cos^ a)y = 0 


( 9 ) 


da 


and 


- (bej. - Sj, cosh^ O y = 0 


( 10 ) 



Similarly the integral in equation (7) would satisfy equations (9) 
and (10) if bej. was replacea by bOj.. Now, since the integrals satisfy 
the separated equations, they must be proportional to products of the 
solutions of these equations. Also, the integrals are periodic in i5, 
so the product solutions must have this periodicity. Consequently, 


I 


2n 


e-inRse (Tj)dTj = pe Ser(a)jer(|) 



-inR 

e SOj.(T])dT) = pOj.SOj.(j^)Jo^(0 


( 11 ) 


( 12 ) 


Jej.(0 and JOp(0 are given as series involving only Bessel functions of 
the first kind, which is \fhy the combinations in (ll) and (12) are 
periodic in i|. 

To evaluate pe^. set a = 0 and develop the left sides of 
equation (ll) in a Bessel function series by using the series definition 
of sej., and the integral definition of the Bessel function (ref. 51 )• 

j (Z) = r ®cos n0 d9 ( 15 ) 

2jti“ ^ 0 

I 

Then evaluate this equation for large values of i by lislng the 
asymptotic formulas for the Bessel functions, equations 4.2.2(27) 
through 4.2.2(50). 

The result is 

pej. = (-i)^ \fSjc (l4) 

To evalmte pOj., first differentiate eqiiation (l2) with respect 


to a. Then proceed as above to obtain the result 
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por = vs; 


(15) 


We may now use ecjuations (ll) and (12) with equations (6) and (?) 
to solve for the expansion coefficients in equation (5). They are 


(1) _ (2) _ Se^.(g) 


nr 


= = 




( 16 ) 


(1) (2) pOj.SOj.(a) 


F'~' = F 
nr ""nr 


(17) 


2N. 


-inR 


Thus the elliptical wave representation of the plane wave is 
= Jej.(5)Sej.(l) + -yl— JOr(S)Sor(Ti)J (l8) 


j 

The expansion for e^“” is simply the conjugate of equation (l8), 
The pressure in the exterior of the cylinder is 



+ ^ D^)Ho^^?]sOrJe^“^ (19) 


Since the external pressiire field is composed of outward -traveling 
scattered waves and the incident plane wave, it may he seen, by 
comparing (19) to k.2.2{j>'^) , that 


= 0, n = 1,2,3 . . 


( 20 ) 
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Therefore 


Pn(5,'l) = Je^(l)jser(l) 

+ Jd^2)ho(2)({) + (.if^ JOr({^SOj.(TiA 

Since the normal velocity is continuous at the cylinder surface, 
equations 4.3(6) and (21) may he used to solve for and 

ci?’ = [inUnr + (-1)''^ Pn Je^(lo)1 

Her Up) r J 


= , *1 

“nr ( 2 )' 

HOr (Sn) 


lnV„- + (-1)' 




The external pressure field is found by substituting (21) through 
(25) into (19) to obtain 

00 00 

P°(5,n,t) = - ^ y |ln[agr(OUnrSer(n) + VnrS°r(’))] 

n=-oc V 



^n ^®nr^®r(’l) "*■ ^®nr^®r(T)j| ® 

(24) 

where 

J 

0 Hey*(e) 



_0 _ r 

%r (2\ 

Her^'Cjp) 

(25) 


0 _ Ho[2)(0 

(2). 

HOr ^ (60) 

(26) 
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^®nr 


l(-i)^^§Tt Sej.(a) 

Nj,He(-)'(5o) 


( 27 ) 


i(-l)^i/BTC So^(a) 

Nr’Ho{.2 )'(So) 


(28) 


For negative values of n, the conjugates of eqxiation (25) through 
( 28 ) are used in the pressure equation. 


4.5 Porous Cylinder Problem 

The interior and exterior solutions have both been given in terms of 
velocity normal to the cylinder surface, that is, in terms of the 
coefficients Uj^ and Vj^. Equation 4.l(ll), vhlch relates pressure 
drop and velocity at the cylinder, will now be used to obtain equations 
which determine these coefficients. 

The pressure drop may be found from the interior solution, 
equation 4.5(9)# and the exterior solution, equation 4.4(24). When 
these are used with 4.l(ll), it becomes 




where 


= Pa 


8q|U 


Of 


n#t) 




8t 


«nr = oSr(5o) 



i 

J®r(®n# * (®n# ^o) 


( 1 ) 

( 2 ) 


"nr 


= - 3 ir( 5 o) = 


JOj. ( Sq, Iq )Ho(^ ) ' ( Sn# So ) 


( 3 ) 
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Fenr = 

and 

FOj^r ■ 


(4) 

(5) 


Since is given in terms of U^^j, and by 4.3(6), equation (l) 

may be \ised to solve for these coefficients. The solution of this 
equation is the fundamental problem for finding the scattering -from a 
porous cylinder. By applying the orthogonality relations for the complex 
exponentials and Mathieu functions, equation (l) may be reduced to 
infinite sets of nonlinear algegraic equations. 
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In equations ( 6 ) and ( 7 ), it must be noted that all Mathieu functions 
and parameters in the equation depend on the parameter Sm 

(eq. 4.2(0)). 

4.5*1 Power Series Equations for Porous Cylinder 
As in the one -dimensional problem, let 


and 


T = R 95 + 6r(qj) 

(1) 

q|~ + . . . 

( 2 a) 

“nr~ + 6u(j> + . . . 

( 2 b) 

V ~ + . . . 

nr nr 

( 2 c) 


then q|^^ is given by 


t) = f y * V^)sOr(l))ei“t 

^ ''LL 


(3) 


n=-oc r=0 v^^“^cosh 25^ - cos 2 t)) 

With these changes, equations 4 . 5 ( 6 ) and 4.5(7) are converted into 

r=0 

r(j) . 


00 


(4) 


) J = 0,1 


+ Fo^f; = (R ^ in^^) ^Nj,(sJV 


(J) 

mr 


r=0 


J 


(5) 
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(6a) 

(6b) 

( 7 -) 

(Tb) 


The power series approximat/ion thus reduces the problem to the 
solution of infinite se^ie of linear algebraic equations. The solution to 
■♦•hese equations are easily approximated by truncating them to obtain 
finite sets of equations. Solving the zeroth-order equations involves a 
matrix inversion and use of the forcing function which depends on the 
incident plane wave. 'The first-order solution uses the same inverse matrix 
and a forcing function which depends on the zeroth-order solution 
(eq. (?a) and (7b)). 

In Appendix A, the Integrals M^j. and Njp are evaluated, and it 
is shown that these integrals ai-e zero unless both Mathieu fvinctions in 
the integral have the same period. Equations (U) and (5) thus divide 
naturally into four sets of equations; one set for each type of 
circumferential Mathieu function. 

* fulfil) “ 

r**0 


(6) 
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00 


N22+l(i®°ta,2Zflui^2Z+l ^®i^2l+l) " ^ 

V Mo 

^”2Z+l,2r+l m,2r+l 
r=0 

(S) 

^2Z+l\,^®^m,22+l^i^L+l + ^°i^2l+l) = ^ 

V No 

/_ "2Z+l,2r+lV2r+l 

(10) 


r=0 


^21+2 21+2 + ^°i^2l+2) = ^ 

00 

^ ^’2Z+2, 2r+2^mf lr+2 

(11) 


r=0 


\ = R + imPA 


(12) 


Equations (8) through (ll) are four infinite sets of linear eqiiations. 
Their solutions may he approximated hy triancating the sets to obtain a 
finite number of equations (ref. 5^)« 

14-.6 Variable Linear Impedance 

There is a remarkable case vf variable material properties vhere 
equation 4.5(1) may be reduced to uncoupled sets of eqxiations. Note 
that if 

Pa = Pa ^2Jt^(cosh - cos 2r\) ( 1 ) 

and 

^ (€Q.|) = R’ V2k2(cosh 2| - cos 2 t)) q|, (2) 

then the procedures described in 4.5 wo\ild result in the solution 

Unr = (5) 

Zn - 

V = (M 

nr ^ ' 

Zjj - inPnr 
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where 

Zi = R' + inp^ (5) 

4.7 Dissipation and Scattering Formulas 
Once the velocity at the cylinder surface is obtained, the pressure 
or velocity at any point may be computed. Some other quantities of 
interest are the incident power, the scattered and dissipated power, and 
the far -field rms pressure. The derivation of these formulas is tedious, 
but routine, and only final results are given here. 

The projected width which the oncoming plane wave "sees” is 

w = h V2(cosh 2| - cos 2 a) , (l) 


and the power which wo\ild pass throia^ this width due to the plane wave 
alone is 


Pj = 2w 


I 


n=l 


The power of the outgoing waves alone is 




00 OC - 

Ul{ 

n=l r=0 ^ 


.(2) 

nr 


Nr(sn) + 


,( 2 ) 

nr 




The rate of energy dissipation within the porous material is 

00 oc 

"" ■ 12 


■I Iff 




U, 


nr 




ntlco 

^®nr^-nrNj.( ®n^ ^°nr^-nr^r(®n) 


'nr 


2 Nl(sn) 


} 


( 2 ) 


(5) 


(M 
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/ 


at large distances from the cylinder where | = » 1, the rms pressure 

is given by 


^00 00 

\n=l ^ ^ r=l 


^nf 2r-2 )^®2r-2( ®n^ T ) + ^^ri( 2r-l )^®2r-l^ ^ ) 

2 ^ 1/2 


(2) (2^ 

■*■ ^®n(2r-l)^°2r-l(®n^^^ “ ^n(2r )^°2r^®n^ 


( 5 ) 



V. RESULTS AND DISCUSSION 


5.1 Linearized Perturbation Theory 

Measured data for the pressure-drop function for a typical sample of 
material are shown on figiire 6 . For comparison, straight lines repre- 
senting acoustic resistances of 0.20, 0.66, and 1.0 are also shown. Each 
linear resistance curve may be tho\aght of as an approximation to the 
material characteristic for a given sound intensity. This intensity is 
related to the value of e. By referring to table I, it can be seen that 
the data of figure 6 are applicable for intensities in the I30 - I60 dB 
range. On figure 2^ these data are plotted on a linear scale for three 
values of e. Since the parameter 5 characterizes the deviation of the 
pressure drop function from the linear curve, it may be seen from this 
figure that 6 may be taken as v/T. 

In other words, for e < 10"^, 6 may be thought of as being large 

in comparison to e, but small compared to 1. This amoimts to saying 
the material nonlinearity becomes important before nonlinearity in the 
wave equation. However, there is no point in adding a 5^ term to the 
power series equations for the material, without also considering 
nonlinearities in the wave equation. Since the equations as derived in 
section 2.1 depended on € being greater than 10"^, there is a definite 
Intensity range where the equations presented herein are valid. In terms 
of the decibel scale given in table I, this range is from about 120 dB 
to 160 dB. The complete theory could be described as having "zeroth order 
wave equations with first order linear boundary conditions in the 
kilohertzian region." This theory is the one used to develop the 
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6? 



Figure 6.- Pressure drop function for typical fiber glass 
plastic porous material. 
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one -dimensional equations In section 3. 5* 5 and the two-dimensional 
equations in section 4.5.1. 

5.2 Nonlinear Perturbation Theory 
If the nonlinear pressure-drop function in 2.2(20) is not expanded 
in terms of the parameter 5, the theory then has zeroth order wave 
equations with nonlinear boundary conditions. Since the boundary 
conditions 2.2(l9), 2.2(20) are accurate to the solutions with this 

theory should be comparable to solutions with the linear perturbation 
theory described in section 5.1. 

The one -dimensional problem depicted in figure 4(b) was used to 
compare the theories and methods of solutions. After obtaining the exact 
solution, the problem was solved numerically by the steepest descent 
procedure of section 3* 3*1. A third solution was obtained by the linear 
perturbation equations of section 3. 3 .3. 

The first harmonic of the velocity at the sample as given by the 
three methods of solution is shown on figure 7* The values of e given 
correspond to incident wave intensities of 137^ 151^ and 157 dB. The 
decrease in velocity amplitude with intensity is due to the general 
steepening of the curves on figure 6. This is called increased acoustic 
resistance. Of course the actual velocity increases with e, since the 
first term in 2.l(l2d) has e as a mxiltiplier. 

The "exact" and perturbation solutions correspond very well for the 
lowest value of e and separate somewhat for larger values. Since the 
error in both solutions is of order this separation is to be expected 

and confirms the earlier statement that the perturbation solution should 



u 



5 X 10”^ 10"^ 


Figure 7 *- Velocity amplitude predicted by exact and 
approximate solutions. 
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_2 

be limited to e < 10 . The correspondence between the steepest descent 

numerical procedure and the other solutions verifies this method, however, 
it seems to be no more accurate than the simpler perturbation procedure. 

For problems other than the special one considered above, the "exact" 
solution is not available and a choice must be me.de between the numerical 
and perturbation method of solution. Since it is far easier to use, the 
perturbation technique seems the logical choice where it is valid. 

5.3 Evaluation of Mathieu Functions 
5 . 3*1 Coefficients and Periodic Mathieu Functions 

Although the coefficients for the periodic Mathieu functions are 
tabulated over the reinge of s used in this paper (ref. ^2), it was 
easier to program the cOTiputer to use the characteristic values to 
compute these coefficients than to read and interpolate the tables of 
coefficients. The method given by Blanch in reference 55 vas used for 
this computation. Results were checked against i:he tables of coefficients 
and were found to check to eight significant figures. With this method 
of computation it is also possible to obtain the coefficients for s > 100 
(which are not tabulated) from the tables of characteristic values for 
large s(s > 100) which are given in reference k2. With the coefficients 
available, the Fourier series in section 4.2.1 could be used for 
computation of the periodic solutions. 

5 . 3*2 Radial Mathieu Functions 

The evaluation of the radial Mathieu functions depends on computing 
the Bessel functions, and Y^. The Bessel f\anctions of the first 
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kind vere computed by Abramowitz' s method, reference 56 . In this method 
a unit value is assigned to Jj^, vhere k is large compared to the 
argument of the Bessel function, is taken as zero and then back- 

ward recurrence is used to find Jn for n < k. The normalization 
condition is then used to adjust the magnitude of the J^s . Yq Is 
obtained from its series representation in terms of J^, Yi is found 
from the Wronsklarv and the rest of the Y^'s are found by forward 
recurrence. Bessel functions computed by this method were accurate to 
nine or ten significant figures for the range of argviments used. 

Since there is little tabulated data available on the radial 
Mathieu f\inctions, these were calculated by two different methods. The 
Bessel function product series given in section 4.2.2 were used and the 
series given by equations 3*05# 3*0*+# 3*15# and 5*16 in reference 42 
were used for comparison. The latter series Involve Bessel functions 
with argument {>Is cosh 0* Each series was checked at the origin (S=0) 
against the Joining factors which are tabulated in reference 42. Also, 
each series was checked against the table of Barakat, Houston, and Levin, 
reference 50 # and against the tables by Blanch and Clemm, reference 57> 
for s = 4. This was the upper limit (on s) of the available tables and 
the lower limit on s for this study. 

The radial functions of the first kind were the most difficult to 
evaluate. Hie product series for JeQ(4,0) checked eigainst the NBS 
tables to seven significant figures. However, this accuracy decreased to 
one signlflceint figure for JCg(4,0), and functions of higher order could 
not be computed by this method. The accuracy of the product series for 
the Je(s,0 and Jo(s,0 functions Increased with s. JeQ( 100,0) was 
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accurate to nine figures, and Je 2 ii^( 100,0) retained a two figure accuracy. 

The Bessel function series proved to be a superior method for com- 
puting Je and Jo. JeQ(^,0) was given to seven figures by this method. 
Although this was not an improvement over the product series it was fo’ond 
that Je5(4,0) was accurate to six significant figui^.*s, \diich was an 
improvement of five significant figures over the product series. Je8(^,0) 
was accurate to five significant figures, but Je 3 ^Q(J'-, 0 ) had only one 
significant figure and JeiI+(^>0) could not be computed. The 

accuracy of the Bessel series for Jej.(s,0) also improved with increasing 
s. JeQ( 100,0) was accurate to six figures, and Je]^^( 100,0) was accurate 
to eight figures. Note that the higher order functions are more accurate 
than the lower order functions when s is large. This txirns out to be 
a very important property, since, in the scattering problem, the low- 
frequency (s=4) solutions require only the lower-order terms for 
convergence. The higher frequency solutions (s=100) require more terns, 
but it is no problem to compute these for the large values of s. 

The radial functions of the second kind were computed by the product 
series. Nej.(4,0) and NOj.(4,0) were accurate to eight significant figiires 
for 0 < r < 15 . The accuracy of the product series for Ner(s,0) and 
NOj,(s,0) decreased with s. For example, Neo(lCX),0) was only accurate 
to three significant figures. NeQ^I^( 100,0) had six figures accuracy; a 
loss of only two places. 

The accuracy of both methods of computation improved with distance 
from the origin. This was established by comparing the results of each 
series and assuming that, if they were the same they were correct. 
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Larger values of 5 correspond to smaller values of cylinder eccentricity 
thus it is generally easier to evaluate the radial functions on a nearly 
circular cylinder than on a nearly flat cylinder. 

In the problem of scattering from a porous cylinder, the product 
series were used to evaluate the functions of the second kind and their 
derivatives. The Bessel series were used to find the functions of the 
first kind, and the derivatives of these functions were found from the 
Wronskian relations 4 . 2 . 2( 58 , 39 ) • 

The magnitude and phase of the complex radial Mathleu functions, 

Hcj. and HOj., could be computed to eight significant figures for s=4, 
and six figures for s=100. This is somewhat surprising in view of the 
above discussion of the accuracy of the functions of the first and second 
kinds which are, respectively, the real and imaginary parts of the complex 
function. Luckily, it turns out that the Inaccurate part of the complex 
number is much less in magnitude than the accurate part, so that the 
complex representation of the number is fairly good throughout the range 
of 5 . 


Evaluation of Coefficients and Integrals 
The evaluation of the coefficients in the scattering eqxiations of 
section 4.5 depends directly on the evaluation of the Mathleu functions. 
The coefficients given by 4.5(2, 3 ) are only as accurate as the derivative 
of the radial function of the first kind. As noted in the previous 
section, this could be computed to five significant figures for s»4 and 
orders up to eight. At s»100, the function of the first kind could be 
determined to six figures for all orders, but the function of the second 



74 


kind could only be obtained to three figures for the lower orders . Since 
the Wronskian was used to solve for the derivati '•e, the derivative is 
limited to three figure accuracy. Improved accuracy for this derivative 
may have been obtained (at s=100) by differentiating the Bessel series for 
the function. Note that the accxiracy problem is associated with the 
waves within the cylinder. The coefficients for the exterior field, 
4 . 4 ( 25 , 26 ), may always be found to six significant figures. 

The nonhoraogeneous terms, 4. 5(4, 5 ) may be computed to six significant 
figures since they involve only the complex radial functions and the 
periodic functions. The normalizing factors; Nj. and are computed 

from the coefficients Dej., DOj. and so are accurate to eight figures. 

Computation of the coupling coefficients, M^j. and is 

discussed in appendix A. The series presented there were used to 
evaluate these. The integrals, 12^, were computed to eight significant 
figures to obtain eight-figure accurewy for the coupling terms. 

The nonhoraogeneous terms in the first-order perturbation solution 
were computed by numerically taking the double Fourier transform of 
r^|(p,t)^ in equations 4.5.l(7a,7b) and expressing the periodic functions 
in ■*'hese eqmtions by their Fourier series so that term-by-term integra- 
tion could be used to find Fe^^^ and accuracy of this was 

checked by solving a linear problem by two methods. The problem \^ere 
^ first solved. In the second method, 

^ = 0-9q{, and r{q^ = O.lq^ was used. Since the two problems 

are physically Identical, the solutions must correspond if the numerical 
work is correctly programmed. It was .'ound that the solutions did match 
to two significant figures, \rtilch is what shovild be expected for 6=0.1. 
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5*^ Scattering from a Porous Cylinder 

5.^.1 Linear Boundary Conditions 
In order to get some general Information about the effects of 
different parameters, solutions were first obtained with linear boundary 
conditions. Focal distances, h, of 2, 4, 6, 8, and 10 were used to 
determine the effect of frequency. These values correspond to cylinder 
lengths of about l/2 wavelength to over 3 wavelengths. Solutions were 
made for acoustic resistances of 1, 2, euad 4. The angle of Incidence was 
varied In 30° Increments from 0° to 90°, and the radial coordinate, 5, 
was taken as 0.2, 0.4, and O. 8 O. Tables II through VI give the ratios 
of scattered and dissipated power to Incident power for these cases. 

Tlie ratio of scattered power to Incident power, SR, varied from 
about 0.04 to over 1.0. It was a mlnlm\im for grazing Incidence and maxi- 
mum for normal Incidence and Increased monotonlcally with resistance. 

As the resistance becomes Infinite, the SR must approach the value for 
the hard cylinder. A computation was made for the case of the hard 
cylinder to check against the work oi Tarakat, reference 28. It was 
noted that the values agreed for normal Incidence but differed by a 
factor of 2 for grazing incidence. Although it has not been possible to 
find the reason for this discrepancy, a possible explanation could be the 
errors which are apparent in table I of Barakat's paper. 

The ratio of the dissipated energy to the incident energy, DR, must 
have a maximum value between R=0, and R=oo. Hie DR depends on the 
product of the pressure drop and velocity at the cylinder surface. For 
zero resistance, the pressure drop Is zero, and for the Infinite 
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TABLE II.- SCATTERED AND DISSIPATED POWER RATIOS, h = 2.0 
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TABLE IV.- SCATTERED AND DISSIPATED PCWER RATIOS, h = 6.0 


R 

1 

2 

4 

1 

a 

SR 

DR 

SR 

DR 

SR 

DR 


0 

0.065 

0.581 

0.161 

0.846 

0.288 

0.984 

0.2 

30 

0.165 

0.670 

0.595 

0.852 

0.732 

0.835 

60 

0.520 

0.758 

0.666 

0.815 

1.072 

0.674 


90 

0.390 

0.747 

0.760 

0.773 

1.165 

0.621 


0 

0.181 

0.590 

0.377 

0.728 

0.599 

0.697 

0.4 

30 

0.236 

0.628 

0.496 

0.754 

0.828 

0.717 

60 

0.358 

0.657 

0.695 

0.685 

1.076 

0.592 


90 

0.376 

0.722 

0.740 

0.742 

1.137 

0.596 


0 

0.272 

0.607 

0.540 

0.738 

0.865 

0.743 

0.8 

30 

0.291 

0.632 

0.584 

0.73** 

0.932 

0.700 

60 

0.340 

0.637 

0.661 

0.713 

1.031 

0.649 


90 

0.377 

0.599 

0.706 

0.661 

1.074 

0.605 
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TABLE V.- SCATTERED AND DISSIPATED POWER RATIOS, h = 8.0 


a 

SR 

DR 

SR 

DR 

SR 

DR 

D 

0.0692 

0.454 

0.171 

0.636 

0.308 

0.710 

30 

C.17^ 

0.663 

0.415 

0.842 

0.773 

0.829 

60 

0.559 

0.667 

0.708 

0.716 

1.109 

0.611 

90 

0.440 

0.609 

0.8l4 

0.596 

1.207 

0.469 

0 

0.122 

0.509 

0.267 

0.704 

0.467 

0.773 

50 

0.261 

0.596 

0.531 

0.715 

0.869 

0.704 

60 

0.543 

0.689 

0.682 

0.748 

1.072 

0.658 

90 

0.424 

0.591 

0.785 

0.585 

1.165 

0.466 

0 

0.188 

0.546 

0.411 

0.700 

0.726 

0.738 

50 

0.301 

0.596 

0.59^ 

0.692 

0.945 

0.671 

60 

0.337 

0.637 

0.660 

0.713 

1.033 

0.648 

90 1 

0.540 

0.658 

0.668 

0.724 

1.042 

0.657 
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TABLE VI.- SCATTERED AND DISSIPATED POWER RATIOS- h = 10.0 


a 

SR 

DR 

SR 

DR 

SR 

DR 

0 

0.085 

0.554 

0.200 

0.807 

0.554 

0.978 

50 

0.186 

0.631 

0.432 

0.795 

0.790 

0.795 

60 

0.384 

0.595 

0.757 

0.615 

1.131 

0.520 

90 

0.429 

0.648 

0.800 

0.664 

1.200 

0.548 

D 

0.193 

0.582 

0.414 

0.735 

0.688 

0.756 

30 

0.256 

0.587 

0.553 

0.694 

0.885 

0.656 

60 

0.37^+ 

0.601 

0.721 

0.630 

1.106 

0.530 

90 

0.381 

0.717 

0.748 

0.739 

1.150 

0.596 

0 

0.289 

0.563 

0.578 

0.650 

0.920 

0.654 

30 

0.304 

0.601 

0.608 

0.686 

0.971 

0.648 

60 

0.31*5 

0.635 

0.680 

0.689 

1.064 

0.608 

90 

0.338 

0.716 

0.687 

0.766 

1.087 

0.647 
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resistance the velocity is zero, while the pressure differential must 
remain finite. It may be seen from the tables that the optimum resistance 
is usually around 2, except for grazing incidence on a thin cylinder where 
it is near 4. This large value of optimum resistance is a surprising 
result, since in one dimensional problems, a resistive impedance of 1 
gives the maximum absorption. For thin cylinders, the DR was an 
increasing function of the angle of incidence for R=l, and a decreasing 
function of the angle for R=4. Another Interesting result is that the 
maximum value for the DR is Just slightly larger than one for grazing 
incidence on a thin cylinder. Thus, it may be concluded that the maximum 
percentage of energy which can be removed from sound propagating in a 
duct is equal to the percentage of duct area blocked by the splitter 
rings and spokes. 

In the case of normal incidence, the dissipation by an elliptic 
cylinder should be roughly comparable to the dissipation by infinite 
parallel planes of porous material. If the length of the minor axis is 
held constant while the major axis is made large, then two sides of the 
cylinder will become nearly straight and parallel. This comparison is 
illustrated by figure 8. The data on this figure are for 5=0.2 which 
corresponds to a cylinder whose major axis is roughly five times the 
minor axis. Dissipation ratios for the use of normal incidence are 
plotted against h, which may be Interpreted as a frequency-parameter. 

The cvirves on figure 8 give the ratio of dissipated to incident energy for 
parallel planes whode spacing is equal to the cylinder's minor axis 
length. The maximum and minimum dissipations occur at very nearly the 
same frequencies for both the cylinder and parallel planes, and the 



R = 

□ R = 2 Elliptic cylinde 

O R = 4 

R = 1 



h 


Figure 8.- Comparison of dissipation by elliptic cylinder 
and ijarallel planes. 
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maxima have nearly the same values. Data for grazing incidence are also 
shown on figure 8, however, instead of the DR, the ratio of energy 
dissipated to the energy in the normally incident wave is shown. This 
emphasizes the fact that only about one-fifth as much energy is 
dissipated from the grazing wave as from the normally incident wave. 

Again the maxima and minima occur at frequencies which correspond to the 
normal incidence case, \rtiich indicates that the minor axis length is the 
critical dimension for grazing wave energy dissipation. The maximxmi 
dissipation occurs when the minor axis is about l/k wavelength. 

'' ^,k.2 Nonlinear Boundary Conditions 

It has been shown in reference 5 that high Intensity waves tend 
toward a sawtooth form at large distances from the source. Because of 
this, the sawtooth has been used to illustrate high-intensity waves 
scattered by an elliptic cylinder. A five-term series approximation to 
the sawtooth wave was used, which gave five time harmonics of the 
fundamental frequency. Figure 9 shows the far-fleld rms pressure 
patterns for the linear solution and the first order power series 
solution for a grazing wave. 

The value of e was 10“^, which corresponds to a sound Intensity of 
about l60 dB. With this intensity, it can be seen from figure 2 that an 
acoustic resistemce of 1 should be used for the linear solution. 

Figure 9 shows that the principle effect of the nonlinearity has been to 
increase the scattered pressure. This change in scattering should be on 
the order of 10 percent, since 6 is 0.1. The curves on figure 9 are 
consistent with these conditions. The largest difference occurs near 




Figure 9*- Grazing incidence 
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the sideline, with an increase of about 50 percent at tj= 85°. Figiire 10 
shows the amplitude of the velocity as a f\mction of position on the 
cylinder surface for the grazing incidence case. The peak velocities 
occxir at the ends of the cylinder, where their value is about 0.95* This 
shows that the choice of e and the resistance was realistic since the 
velocity varies between 0 and 1. 

^so shown on figure 10 is the maximum deviation of the actual 
pressure drop from the pressxire drop predicted by the linear resistance. 
Since it is always less than 0.2, it may be regarded as a perturbation. 
Again, this Justify' s the use of the power series technique for this 
problem. 

Far-fleld pressure patterns for the normal incidence case are shown 
on figure 11. Although the only difference between this case and the 
previous one is angle of incidence, the large deviation between the 
patterns suggests n violation of the assumptions necessary for power series 
analysis. An inspection of the velocity amplitude at the surface for this 
case reveals that it has been underestimated by an order of magnitude. 
Consequently this solution is not valid. A better estimate of the true 
nonlinear solution would require a value of € of 10“^, but this is 
beyond the applicable range of this method, as was pointed out earlier. 
Thus, a solution for a l60 dB wave may be found for grazing incidence, 
but solutions must be restricted to lower intensities for normal 


incidence. 




8 



VI. CONCLUDING REMARKS 


Integral conservation theorems have been used to derive a set of laws 
which govern the discontinuities in mass, momentum, and energy flows at 
a thin sheet of porous material. Each of these laws involve two unknown 
functions ■vrtiich must be evaluated experimentally. The simplest reasonable 
approximation to these general laws is made by ass\milng continuity of 
normal velocity at the porous surface, and a pressure discontinuity which 
is a nonlinear function of only the normal velocity. Since this pressure 
drop is a function of normal velocity, it may be found from a simple 
steady-state test. It is useful to approximate this nonlinear function 
by a linear function \rtiose slope depends on the intensity. Deviations 
from this c\irve may be represented by a parameter which also 
depends on the intensity. This approach has been shown to be useful for 
intensities up to about l60 dB. 

The spectral analysis formulation of a one -dimensional problem 
results in infinite sets of nonlinear eqixations. These may be sol/ed by 
an optimization procedure, such as the steepest descent method, for 
arbitrary nonlinearity or by power series techniques for small nonlinearity. 
The power series technique is also useful when it is combined with a 
quasi -linear approximation to the actual nonlinear function. 

The temporal formulation of the one -dimensional problem results in 
a nonlinear integral equation. Numerical approximations to this Integral 
equation again give finite sets of nonlinear algebraic equations to be 
solved. 
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Two exact solutions of the one -dimensional problem may be found. 

These solutions are for a porous sheet mounted in an infinite tube and 
for a sheet with a rigid quarter-wavelength backing. Both solutions are 
given in terms of the inverse of a nonlinear algebraic function. 

The exact solution for scattering from a porous elliptic cylindrical 
shell with linear acoustic resistance may be found in terms of Mathieu 
function series. Coefficients of these series must be foimd by inversion 
of four infinite matrices. Elements of these matrices depend on integrals 
of the periodic Mathieu functions and on the radial Mathieu functions. 
Series expressions for the Integrals have been found. The Mathieu 
functions may be evaluated in terms of known trigonometric and Bessel 
function series. By use of highly accurate modem computers, it has been 
possible to evaluate these functions for a range of parameters where no 
tables are available. Thus, it has been possible to present results for 
a range of frequencies where only formal solutions exist in previous 
literature. 

For a special case where the acoustic resistance is a certain 
function of position on the cylinder, the infinite matrices reduce to 
diagonal forms, so that simple closed-form expressions in the series 
solutions are available. 

It has been possible to obtain numerical results inside the cylinder, 
and outside of the cylinder in both the near -field and the far -field. 

Near -field solutions have been used to obtain the power dissipated in the 
porous shell. 

The power series method has been used to obtain approximate solutions 
for scattering from a cylinder with nonlinear acoustic resistance. It 
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has been shown that there may be significant differences between the linear 
and nonlinear solutions at high intensities. 

The exact solution for scattering from a cylinder with large non- 
linearity depends on solving doubly-inf inite sets of nonlinear algebraic 
equations. It has been fo\md that ccmputer-time for a power series 
solution is roughly 100 times as large as for a linear solution. Because 
of this, it is impractical to solve the problem of large nonlinearity in 
the spectral domain \mtil more efficient computational programs can be 
developed. 
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DC. APPENDDC 


Evaluation of Coupling Coefficients 


The integrals defined in equations 4.U(8) and 4.4(9) niay be evaluated 
by expanding the radical in the integrand in a series and then integrating 
term by term. The integrals to be considered are 


M 


It 


°hfr 


2it Se^(Ti)Sej,(Ti)dTi 

0 \|2(cosh 2|- - cos 2 t]) 


and 


= f / 


■2n Soi{T\)So^{r\)dr\ 

0 \|2(cosh 2 |q - cos 2 t)) 


( 1 ) 


( 2 ) 


The radical in the integrand of equations (l) and (2) may be 
represented as a series of Gegenbaur polynomials. Rainville (51) gives 
the generating function definition of these polynomials as 


(1 - 2xt + 


00 



n=0 


The denominator of equations (l) and (2) may be written as 


(5) 


|2(cosh 2^0 - cos 2 = ^e^^°(l - 2 cos 2 t| e”^^° + (4) 


Therefore 


[2(cosh 25 q - cos 2 tH ^ = e "^o Y C^^^(cos 2 t|) e 


(5) 
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The Gegenbauer polynomial in equation (5) may be given as (ref. ^1, 
p. 285). 


cV2(cos 2ti) = y 

” ^ kl (n - k)J 


k=0 


( 6 ) 


Equations (5) and (6) combined give the desired expansion for the 
radical 


[2(cosh 2 ^q - cos 2q)J 


1/2 -% 


I 


= e ” ^ ® ° 

n=0 lt=0 


r (l)k©n-k 


k! (n - k): 


(7) 


From equation (7)^ it may be seen that the integral 


cos mT] dn 

^ra = / - - ■ 

0 \/2(cosh 2|q - cos 2tj) 


( 8 ) 


is zero for all odd values of m, consequently, the only integi'als which 
it is necessary to evaluate are of the form 


^2rri 


r 2jt 

~^Q \ 


cos 2mTi dT) 


\f2(cosh 25o - cos 2 t)) 


(9) 


Substituting equation (7) into equation (9) gives 


l2m = e 


-if 


>y r%os a(n - 2k). 

Lj Z- kl (n - k)! Jq 

n=m kaO 


cos 2mii dT) ( 10 ) 


kaO 

k! (m + 2p - k) I 0 


l2m*e (?)^(^ ( l atSP-h) f cos 2(m + 2p-2k)T) cos 2mi) dr) 

p=0 k* 


(n) 
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In equation (ll), nonzero terms occur only when k=pork=p+m, 
therefore, 


^2m ® 


p=0 


-(2m+l)^o^ \2/p\2/p+m 
L p! (p + m)! 


( 12 : 


Since the factorial terms in equation (12) form a decreasing sequence, 
the error due to truncating the series is 


m fi\ -i*(N-i)eo 


N! (N + m)! i+lo 

A relative measure of the error is obtained by taking the ratio 
of ( 15 ) to the first term in the series for Ip^. 


( 13 ) 


E T < 
rel 


r (2m + l)(2m + 3)- » -(2m + 2N - if l e 
L(2m + 2)(2m + i+). . .(2m + 2Nl) 


-i^(N-l)r 


W 


Since the first term in the inequality (l4) is always less than 1, 

-45o(N-l) 

^ei rr 




Then if 


and 


^rel 


= 10 


.-P 


-R 


( 15 ) 

( 16 ) 
( 17 ) 


an estimate of the number of required terms is 


N > 10^(P > R) 


( 18 ) 
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Thus, for eight significant figures (P = 8), 90 terms are required 
for Iq = 0.1, and 1000 terms are required for §o = 0.01. With N 
determined hy (l8), we put 


N 


l2m = 2,t e 


-(2m+l)lo^ 


(B.® 


p+m 


-^P5( 


L p! 

p=0 


(p + ni)l 


(19) 


Once the integrals are evaluated, they may be used to evaluate 

the integrals in equations (l) and (2). First, consider the product of 
the Mathieu functions in the integrand of equation (l). If sOj^ and 
sej. both have period n, their product may be written as 


^®2L^®2r 


■I !■ 

k=0 j=0 


(2L) _(2r) 

2J 


2k 


De 


cos 2kq cos 2jq 


(20) 


cos 2kq cos 2jT) = i cos 2(k - j)''! + ^ cos 2(k + J)t] 


Note that each term in the series is an even m\;iltlple of t). 
Consequently, 


^L,2r 


OC CO 



k=0 k=0 


2J \ 2 


( 21 ) 


The product of Mathieu functions with different periods gives a 
series whose terms are all odd multiples of t). It was noted in 
equation (8) that the integral of the^e terms is zero, therefore. 


^21/4-1, 2r = 0 (25a) 

^2L,2r+l = 0 


(23b) 
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The product of even Mathieu functions with period 2n is 
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Therefore 
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A similar evaluation for integrals involving odd Mathieu 


gives 
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functions 
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